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Abstract 

We investigate the quantum dynamics of systems involving small numbers of strongly interacting 
photons. Specifically we develop an efficient method to investigate such systems when they are 
externally driven with a coherent field. Furthermore, we show how to quantify the many-body 
quantum state of light via correlation functions. Finally we apply this method to two strongly 
interacting cases: the Bose-Hubbard and fractional quantum Hall models, and discuss an imple- 
mentation of these ideas in atom-photon system. 
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Strongly interacting photonic systems provide new avenues for examining quantum sim- 
ulation, topological quantum computation, and many-body states of matter (lH5|. Inspired 
by analogies to electronic systems, these studies focus on ground state properties. How- 
ever, a photonic system is naturally an open driven system. Therefore, the most relevant 
approach to understand and manipulate many-photon states involves understanding the 
non-equilibrium dynamics in such systems 6l4l l| . In particular, this difference raises key 
questions: given the presence of photon loss, how does one prepare a photonic state with 
many-body features? What is the manifestation of important properties such as incompress- 
ibility and collective effects, when the system is coherently driven with a laser field rather 
than coupled to a thermal bath? 

In this article, we address these questions by studying a driven system of strongly in- 
teracting photons and evaluating physical observables that display quantum many-body 
signatures. We focus on a two-dimensional lattice of interacting photons with an effective 



gauge field 



12 



15] . In the presence of strong interaction (nonlinearity) on each site, the 



system maps into the bosonic fractional quantum Hall 



16Ml8| model. Such nonlinearities 



have been experimentally shown at the single site (resonator), both for optical 19M2l| and 



microwave 



22 



23| photons. We demonstrate that by weakly driving the system, a few 



photon Laughlin state can be prepared. We introduce experimentally-relevant observables 
such as the correlation function of the common-mode (the common-mode) to investigate the 
response of the system. Furthermore, we present a scheme to adibatically prepare such state 
using many-photon Fock state and compare it to a driven scheme. 

The key idea underlying our approach to the driven scheme is to generalize the theo- 
retical technique of weakly driven cavity-QED system (involving one atom interacting with 
one photon) to the many sites and many-photon regime. We follow Carmichael et al. 24 1 
who showed that when an optical resonator with strong nonlinearity is weakly driven, one 
can truncate the Hilbert space up to two-excitation states, reduce the exact master equa- 
tion description to an effective Schrodinger equation description. Most importantly, the 
"quantum jumps" do not contribute in the correlation functions. In particular, while the 
one-photon state is intact under nonlinearity, the two-photon component exhibits bunching 
or anti-bunching features. Similarly, in a system of many sites, the photonic state at each 
photon-number manifold reorganize themselves according to the interaction. In a dilute 
lattice with Nj, magnetic flux quanta and strong interaction for a fixed number of bosonic 
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particles (N p h), the system is expected to have fractional quantum Hall states (Laughlin 



type) at filling factors v = Nph/N^ = 1/2, 1/4, ... 25|, |26|. We demonstrate that when an 
optical system is driven with a weak coherent field, which has Poissonian distribution of 
photon number, the system forms Lauhglin state in a photon-number (N p h) manifold which 
corresponds to the bosonic Laughlin filling factors, at specific pump frequencies. We show 
that measuring the A^-body correlation function reveals the existence of such state. Fur- 
thermore, we present an alternative adiabatic method to prepare such a state for larger 
photon number and compare the two methods. 

We consider a 2D interacting photonic system which has the Hamiltonian {h— 1): 

H - — T \ n f, i27Ta y i p,\ ft , p-«27roy 

-"sys ° / j u, x+l,y u, ^,y° ~ x,y u, x+l,y c - 

+ & x ,y+l® J x,y ~\~ ^x,y^x,y+l 4" -^frcc + H{ ni . (1) 



where a\ y is the creation operator at site (x,y), J is tunneling rate between resonators, a 
the effective magnetic flux per plaquette (total magnetic flux is = aN x N y ), and Hf ree = 
Y2x y u o a x y a x,y We take an on-site interaction term of the Kerr-type: H mt = U a\<ii(a\ai — 1) 
where the index stands for the site i = (x,y). In the absence of the magnetic field (a = 
0), the Hamiltonian describes the Bose- Hubbard model. We return to the discussion of 
implementation of such Hamiltonian extending the scheme proposed in Ref. [13] later in the 
article. 

To include loss and driving, we use the stochastic wave-function approach 



24 
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The coherent drive is applied uniformly; its effects and that of the associated loss can be 
described by the non-Hermitian term H pump = ^ i K(5{e lu)pt a i + e~ lulpt a\) — iKa\a,^ where k 
is the coupling rate to the resonators, f3 is the amplitude and u p the frequency of the drive 
field. In the rotating frame of the pump field, the effective Hamiltonian of the driven system 
is: 

H cS = H sys + nf3 Y^ifli + at) - (A + in) ^ n t (2) 

i i 

where the pump detuning A = u p — ojq takes the form of a chemical potential. Since the 
system is open, in the absence of the pump ((3 = 0), the system will be in the vacuum state. 

We generalize the quantum-jump picture for evaluating the correlation functions j^] to 
many-photons and many-modes. The evolution of the system is governed by the effective 
Hamiltonian (Eqj2]) and the corresponding quantum jump operators (<3j). In particular, in 
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the weakly excited system (,8 < 1), the metastable state of the system can be perturbatively 
written as: 

|tf> ~ |0) + 0{$)\1) + G(/3 2 )|2) + ...0(P n )\n) + ... (3) 

where \n) = ^2 Cj 1 ...j n aj l ...a|jO) represents a state in the n-photon manifold of the lattice 
system. This state is the eigenstate of H e ff with the smallest imaginary eigenvalues, i.e. it is 
mostly the vacuum state. All other states have at least one photon, and therefore, they decay 
rapidly into this state. When a photon decays from any site, the system undergoes a quantum 
jump. These jumps occur at a rate kO{[5 2 ), and the system takes a state of the form: 
|*;) =a i |*)/(|o i |*)|) ~ \0) + O{f3)\l') + O{f3 2 )\2') + ...O(f3 n )\n') + .... Similarly, the system 
can undergo a two-photon jump with a slower rate nO(f3 4 ). Since the system is continuously 
pumped, it is restored back into the steady state with a relatively fast rate (k) . Therefore, the 
density matrix of the system can be formally written as: p = \ty)(fy\ + 0(/3 2 )pi + 0(/3 4 )p2 + ..., 
where pj stands for the density matrix after "j" consecutive jumps. For a single jump we have, 
Pi = Now > we evaluate the n-body correlation 

function of an arbitrary operator d which is a linear superposition of the site operators (ctj). 
In particular, we are interested in = {pd) n d n ). Using the above picture, this correlation 
function can be perturbatively written in powers of pump amplitude: 

G (n) = 0{f3 2n )(n\d r d n \n) + 0(p 2n+2 )(n'\£ n d n \n') + ... (4) 

Therefore, if we are interested in the n-photon manifold, the metastable state \^f) is sufficient 
for evaluation of any n-body correlation function and the corrections due to quantum jumps 
can be ignored. 

In particular, for a two-particle case, we define the two-body observables to characterize 
the deviation from the classical regime. For a single resonator this deviation is characterized 
by the equal time second-order correlation function as g^> = (a] 2 a 2 ) / (d\di) 2 . This quantity 
is useful in characterization of cavity QED experiments. However, this observable can not 
encapsulate the collective effects in the system. In particular, in the presence of strong 
interaction (U ^> J), such a quantity is always less than one regardless of the collective 
features of the entire system. Instead, we consider a collective observable which is the 
second-order correlation function of the common-mode = ^= 9cm = ^tg^ • ^ ms 

observable is particularly interesting since we are exciting all the resonators the same way, 
and therefore, this mode is primarily excited. Such quantity can be obtained by measuring 

4 



the 2nd-order correlation function of the far-field light emitted from all the resonators. 

To numerically find we consider a truncated Hilbert space corresponding to at most a 
few particles and find the eigenstate of H e g with the smallest imaginary part of its eigenvalue. 
Note that in contrast to grand canonical ensemble - where we minimize (H — fiN) - here 
we find the steady state of the system as a function of the pump field. Such approach 
allows us to consider larger lattices which are otherwise inaccessible with the density matrix 
approaches [17] . 

Using the technique described above, we study the driven system of interacting photons 
with the Hamiltonian of Eq.flTJ. First, we consider the case of hard-core bosons (U ^> J), 
and investigate the response of the system as a function of the pump field frequency (A). 
For simplicity, we only consider the case of v — 1/2. The input field consists of a Poisson 
distribution of photons. When photons are injected at the frequency corresponding to the 
Laughlin state at the A^-photon manifold, photons reconfigure themselves and form a wave 
function which corresponds to the Laughlin state. The remarkable overlap of this photonic 
state with the Laughlin wave function in the A^^-photon manifold is shown in Fig. 1(a). 
Note the frequency required to be resonant with the Laughlin state is at the vicinity of the 
free photon state (Hofstadter's spectrum). In the limit of large system (N x N y — > oo), and 
dilute magnetic field (a C 1), these two frequencies coincide since the Laughlin state is the 
exact ground state of the Hamiltonian in the continuum limit. For numerical simulations, 
we have used the discrete version of the Laughlin wave function on the lattice with torus 
boundary condition [26]. 

Around the resonance, we observe the suppression of the correlation function of the 
common-mode. The reason behind this suppersion is that the external pump is coupled 
differently to the single particle manifold and Ap^-photon manifold, corresponding to the 
Laughlin filing factor. We note that the energy of the single particle state and the Laughlin 
state per particle is exactly equal to each other in the continuum limit, and the previously re- 



ported discrepancy is due to the finite size effect 



151 ] . The direct experimental verification of 



the Laughlin overlap is a difficult task which requires number post-selection {N p h) and state 

( N X Ny\ 

tomography in a Hilbert space with dimension . However, the common-mode 

\N ph J 

correlation function can be obtained by using conventional quantum optics measurements. 
Now, we relax the hard-core constraint and investigate the same observables. In the weak 
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Figure 1: Overlap with the Laughlin wave function (y = 1/2), and the correlation function of the 

(2) 

zero mode {g CM ) are shown as a function of: (a) the pump frequency for hard-core bosons (b) the 
interaction strength for A = — 3.36J. The simulation is for a 6x6 lattice, torus boundary condition. 
K = .01J, (3 = 0.01 and Na, = 4, and we have evaluated the overlap with the Laughlin function 
in N p h = 2 manifold. (c,d) are similar to (a,b) for N p h = 3, = 6, A = — 3.095J and the 
corresponding correlation function {g CM )- All calculated quantities are dimensionless. 



interaction limit, the system approaches the classical response, as shown in Fig. 1(b). In 
the absence of interaction, using transport measurements - varying the pump frequency and 
measuring reflection/transmission- one recovers the Hofstadter's butterfly spectrum 13], 
but regardless of the pump frequency, the correlation function remains equal to one. Similar 
behavior was observed for N p h = 3, as shown in Fig.l (c,d). 

In other to clarify the connection between zero mode correlation and the collective nature 
of the system response, we investigate the driven photonic Bose- Hubbard model (Fig. 2). 
In the limit of weak interaction, the system behaves classically and the correlation function 

(2) 

approaches that of a coherent state, i.e. g CM = 1, as shown in Fig. 2(b). However, 
in the strong interaction limit (U ^> J) the system exhibits significant deviation from a 



6 



1Cf 



10 



10" 



10 



10" 



(a) 

i 
i 

~ ~ v ' 

V 1 

* i 
» i 


'i 

1 1 
\ 

\ 

\ 


---g (2) 


•i 
>i 
i 





10 



10" 



10" 



10" 



(b) 



-4 -3.9 -3.! 

(0 /J 
P 



-- g 



(2) 
CM 



.01 1 100 
U/J 



Figure 2: The correlation function of the zero mode {g CM ) are shown as a function of: (a) the 
pump frequency for hard-core bosons (b) the interaction strength for A = — 4.0J (where g CM is 
maximum). The simulation is for a 6x6 lattice, torus boundary condition and the maximum number 
of photon is 3. n = 0.002 J, (3 = 0.01. 

classical state. This deviation is due to the finite size of the system and can be understood 
in the following way: the system is weakly driven and manifolds with large number of 
photons are weakly populated. Therefore, the effective filling factor (n tot ) / (N x N y ) is small 
and in the presence of a non-zero interaction, one expects the system to be in a superfluid 
regime. However, due to finite size of the system, the common-mode is not completely 
harmonic and the two-photon resonance is slightly shifted. This leads to a deviation of the 
correlation function from unity; using the single-mode approximation, we get an estimate 
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(2) 



We numerically verify such statement by evaluating the correlation function (gc M ) as 
a function of the system size. In the Bose-Hubbard model, as the system size increases, 

(2) 

the correlation function g CM approach the classical limit, i.e. unity, as shown in Fig. 3(a), 
while the correlation function of individual sites is equal to zero. The green curve shows the 
numerical estimate based on the nonlinearity between one- and two-photon manifold lowest 

(2) 

energies, which diminishes as the system size increases. In contrast, in the FQH model, g CM 
remains constant as the system size changes, as shown in Fig. 3(b). Note that the overlap 
with the Laughlin wave function is also constant and remains close to unity. We have 
also performed numerical simulation for two-point correlation function g(i,j) = (alalajdi) 
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Figure 3: The correlation function of the zero mode {9zm) are shown as a function of the system 
area (N x x N y ) for: (a) Bose-Hubbard and (b) Fractional Quantum Hall models. N p hi = 2, K = 
0.04J, (3 = 0.01 and the system is truncated at three photons. 

projected into the iV^-photon manifold, and the results agrees with two-point correlation of 
the Laughlin state. Note that in the general case of n-photon FQH state, one should measure 
n-body correlation function G^ 1 ' = (pcft n d n ), as introduced earlier. Such correlation function 



29j : the photonic mode d is 



can be measured using a modified Hanbury Brown- Twiss setup 
collected, the light passes through n beam splitters and then the state is detected using n 
photodetectors. 

Now we discuss the implementation of the Hamiltonian in Eq.([T]) and the conditions 
to observe fractional quantum Hall states of photons. Recently, there have been several 
proposals to implement the artificial magnetic fields for photons 



achieve strong interaction in coupled resonators systems 
in Ref. 



US- 



12 



and various means to 



Here, we focus on the proposal 



13] which does not require time-reversal symmetry breaking for the implementation 
of the magnetic field. Strong photon-photon interaction -which can lead to photon blockade 
- can be mediated by coupling emitters (e.g., atoms 30|, quantum-dots 3lfl, Rydberg states 



32 



34] for optical photons and Josephson junctions for microwave photons [35]) to the 
resonators. 

Besides the driven method to reach fractional quantum Hall state that we discussed above, 
one can also prepare a Laughlin state by adiabatically melting a Mott-insulator of photons, 
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similar to the atomic method discussed in Ref. |25j, as described in Fig. HI However, this 
requires both preparation of N p h Fock states and photon lifetimes long enough to allow for 
the melting to be adiabatic, making the coherent drive approach preferable. Note that one 
might be able to use the nonlinearity of the system itself to prepare the N p h Fock states of 
photons [36]. 

Regardless of the preparation method, coupling atoms to the photonic system introduces 
loss which can be reduced by detuning the cavity resonance from the emitter transitions 
(A, A' 3> T). As an example case, one can use an ensemble of N- level atoms to mediate onsite 
two-body interaction of the Kerr-type (FigJH(b)) js], which still preserves the propagation 
direction (clockwise or counterclockwise) used in Ref. 13] . In this approach, the optical 
cavity and ensemble enter into a slow-light regime, where the excitations are dark state 
polaritons 



37j | ^x,y oc fia Z)J/ — gyNS Xt y, where Q is the pump field, g is the vacuum Rabi 



coupling, N is the number of ensemble atoms, and S x>y is the spin-wave operator describing 
coherence between two atomic states \a) and |c) (from Fig. H^b)). These bosonic excitations 
lead to an overall increase of dynamical timescales by r] = c/v g ^> 1, the ratio between the 
speed of light and group velocity for the dark state polariton, but they can also interact via 
a self-Kerr interaction with state \d) 38| . For observing a Laughlin state and having a finite 
gap, the effective interaction between photons (U ~ g 2 /A') should be at least comparable to 



26] . These conditions can be satisfied for systems with a large Purcell 



the tunneling rate J 

factor (g 2 /nT > 1). The same criterion applies to implementation of such scheme in the 
microwave domain. 

In conclusion, we have shown that driven strongly interacting photons exhibits interesting 
many-body behaviors and FQH state of photons and their incompressibility can be probed 
by using conventional optical measurement techniques. Investigation of other many-body 
signatures of these states such as their topological properties and fractional statistics and 
preparation of photonic many-body state with reservoir engineering (39)] can be the subject 
of further research. 

This research was supported by the U.S. Army Research Office MURI award 
W911NF0910406, NSF through the Physics Frontier Center at the Joint Quantum Institute, 
CUA, Packard, Darpa and AFOSR MURI. We thank E. Dernier and I. Carusotto for fruitful 
discussions and E. Goldschmidt and S. Polyakov for critical reading of the manuscript. 
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Figure 4: Adiabatic preparation photonic Laughlin states: (a) Atomic ensembles are coupled 
to resonators to mediate interaction. A control field couples internal levels of the atom, shown in (b), 

n 

and provides on-site interaction for photons [3j • (c) Overlap of the two lowest states with Laughlin 
wave function (energy levels relative to the ground state) are represented by dashed (solid) lines, 
respectively. The procedure to make a Laughlin state: (i) Create N p photons in the whole system 
(e.g., by using lambda systems inside the resonator), at this stage a is set to be zero, (ii) Make 
a N' x x Ny superlattice potential V (e.g., by detuning selected resonators) such that the ground 
state gets to the first Mott insulator (N p h = N' x N' y ). (iii) Turn on a single-site potential V per t by 
detuning a cavity (in this case (x,y)=(3,3)). (iv) Turn on the magnetic field to the desired value 
a = Nph/ (vN x N y ). (v) Melt the Mott insulator by lowering the superlattice potential strength 
to zero, (vi) Lower the single-site potential. Three snapshots of lattice potential are shown at 
(iv) and the end of (v) and (vi) steps, respectively. The impurity potential splits the ground state 



degeneracy of the Laughlin state on the torus boundary condition 26] and prevents level crossing 
and sharp changes in the overlap. 
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